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ABSTRACT 

The  evolution  of  forces  acting  on  horizontal  cylinders  subjected  to 
impact  by  an  harmonically  oscillating  free  surface  has  been  investigated 
both  theoretically  and  experimentally. 

The  experiments  were  carried  out  in  a  large  U-shaped  water  tunnel, 
with  cylinder  diameters  ranging  from  3.0  to  6.5  inches.   The  results 
have  been  expressed  in  terms  of  two  force  coefficients.   The  first  is 
the  slamming  coefficient  which  expresses  the  normalized  force  acting  on 
the  cylinder  at  the  time  of  impact;  and  the  second  is  the  maximum  drag 
coefficient  which  occurs  when  the  cylinder  is  immersed  approximately 
1.75  diameters  in  water.   The  slamming  coefficient  was  found  to  be  equal 
to  it,  as  predicted  theoretically.   It  was  also  found  that  the  slamming 
coefficient  may  be  amplified  to  a  value  as  high  as  2tt  through  the 
dynamic  response  of  the  cylinder.   The  maximum  of  the  normalized  force 
in  the  drag  dominated  region  was  found  to  be  approximately  equal  to  2.0. 
It  decreased  with  increasing  Froude  numbers  to  a  value  of  about  unity. 

It  is  recommended  that  further  experiments  be  conducted  with  sand- 
roughened  cylinders  and  the  dependence  of  the  rate  of  rise  of  the  impact 
force  on  the  flow  parameters  be  investigated. 
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I  .   I NTRQDUCT I  ON 

information  concerning  the  forces  acting  on  bluff  bodies  subjected 
to  wave  slamming  is  of  significant  importance  in  Mechanical  Engineering, 
Ocean  Engineering,  and  Naval  Architecture.   The  design  of  structures 
which  must  survive  in  a  wave  environment  is  dependent  on  a  knowledge  of 
the  forces  which  occur  at  impact,  as  well  as  on  the  dynamic  response  of 
the  system.   Two  typical  examples  include  the  structural  members  of  off- 
shore drilling  platforms  at  the  splash  zone  and  the  often-encountered 
slamming  motion  of  ships. 

The  general  problem  of  hydrodynamic  impact  has  been  studied  exten- 
sively EG  motivated  in  part,  no  doubt,  by  its  importance  in  ordnance 
and  missile  technology.   Extensive  mathematical  models  have  been  developed 
for  cases  of  simple  geometry  such  as  spheres  and  wedges,  and  these  models 
have  been  well  supported  by  experiment.   Unfortunately,  The  special  case 
of  wave  impact  has  not  been  studied  as  extensively.   Kaplan  and  Silbert 
[2~]   developed  a  solution  for  the  forces  acting  on  a  cylinder  from  the  in- 
stant of  impact  to  full  immersion.   Da  I  ton  and  Nash  Q3H  conducted  slam- 
ming experiments  with  a  0.5  inch  diameter  cylinder  with  small  amplitude 
waves  generated  in  a  laboratory  tank.   Their  data  exhibited  large  scatter 
and  showed  no  particular  correlation  with  either  the  predictions  of  the 
hydrodynamic  theory  or  identifiable  wave  parameters.   It  appears  that  the 
evaluation  of  the  slamming  effects  with  wavy  flows  is  extremely  difficult 
partly  because  of  the  limited  range  of  wave  amplitudes  that  can  be  achieved 
and  partly  because  of  the  difficulty  of  measuring  the  fluid  velocities 
at  the  instant  of  impact. 
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In  essence,  the  present  research  was  undertaken  to  extend  the  work  of 
Da  I  ton  and  Nash  and  to  improve  the  experimental  method  by  using  a  large 
U-shaped  water  tunnel.   In  addition,  it  was  intended  to  demonstrate  a  cor- 
relation between  theory  and  experiment  by  taking  into  account  the  dynamic 
response  of  the  system.   Specifically,  the  goals  of  the  project  were  as 
fol lows: 

A.  to  examine  the  existing  theoretical  models  in  determining  wave 
slam  forces  on  circular  cylinders; 

B.  to  furnish  data,  obtained  under  controlled  laboratory  conditions, 
about  the  forces  acting  on  smooth  circular  cylinders  subjected  to 
impact  with  an  harmonically  oscillating  water  surface; 

C.  to  determine  the  relative  importance  of  the  inertia  dominated  and 
drag  dominated  forces  during  wave  impact;  and 

D.  to  correlate  such  data  with  respect  to  identifiable  wave  parameters 
such  as  the  Froude  number,  Fr;  the  Keu I egan-Carpenter  number,  K; 
and  the  Reynolds  number,  Re. 
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I  I .   THEORETICAL  ANALYSIS 

The  traditional  approach  used  in  the  design  of  off-shore  structures 
involves  the  classical  Morison  equation  to  determine  the  forces  due  to 
wave  motion.   However,  wave  impact  is  generally  defined  as  the  early 
stages  of  the  penetration  of  a  sol  id  body  into  a  wave  surface.   At  the 
instant  of  contact  the  fluid  in  the  vicinity  of  the  body  undergoes  large 
accelerations  which  give  rise  to  large  forces.  As  the  body  becomes  more 
fully  immersed,  forces  due  to  viscous  drag  and  separation  effects  become 
predominant.   Thus  the  inertia  and  drag  coefficients  used  in  the  Morison 
equation  are  not  constant, and  the  problem  becomes  very  difficult,  even 
for  simple  geometries. 

The  general  case  of  hydrodynamic  impact  is  usually  described  by  using 

incompressible  potential  flow  theory.   For  the  case  of  a  moving  body  with 

mass, M, and  velocity, v,  impacting  a  quiescent  surface,  the  system  momentum 

i s  Mv  .   Neglecting  non-conservative  forces,  the  momentum  of  the  system 

is  unchanged  during  penetration.   However,  the  mass  of  the  system 

increases  due  to  the  fluid  which  is  set  in  motion  in  the  vicinity  of  the 

body.   Also  known  as  "added  mass",  m,  it  results  in  reduction  of  the 

velocity.   Thus  the  system  momentum  after  penetration  is  (M  +  m)  v  =  Mv  . 

o 

The  impact  force  at  any  instant  is  a  function  of  m  and  3m/9t.   The 
sol ut ion, therefore,  requires  knowledge  of  the  added  mass  and  its  time 
derivative. 

It  should  be  noted  that  viscous  effects,  bei ng  non-conservative,  will 
alter  this  somewhat.   Additionally,  high  speed  entry  (as  with  projectiles) 
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may  require  that  compressibility  effects  of  both  the  fluid  and  the  air 
above  it  be  taken  into  account.   However,  for  low  speed  entry,  these 

effects  are  generally  neg lected, and  the  fluid  motion  is  described  by  a 

2 
velocity  potential  <j>  ,  which  satisfies  V  <j>  .=  0.   The  added  mass  term 

can  then  be  calculated.   The  problem  remains  extremely  difficult,  though, 
because  the  free  surface  tends  to  "pile  up"  around  the  body.   This 
results  in  an  unsteady  flow  problem  with  a  time  dependent  boundary  condi- 
tion (the  free  surface).   Approximate  solutions  have  been  obtained  E-O 
for  simple  shapes  (spheres,  cones,  cylinders,  wedges,  flat  plates,  etc.) 
by  assuming  a  constant  plane  free  surface  and  then  applying  a  correction 
for  the  distortion. 

Kaplan  and  Silbert  \JZ~^   developed  the  following  relation  for  the  force 
acting  on  a  horizontal  circular  cylinder  by  a  wave  system  which  propa- 
gates normal  to  it,  (See  Fig.  |.): 


£  =  pg  A.  +  (m  +  p  A. )  r\  +  —  r\2  ( I  ) 

L    3   i  i      8z 


in  which  F  represents  the  force  acting  on  the  cylinder;  L,  the  length 
of  the  cylinder;  p,  the  fluid  density;  g,  the  gravitational  acceleration; 
A.,  the  immersed  area;  m,  the  added  mass  per  unit  cylinder  length;  n.,  the 
instantaneous  height  of  the  wave  surface  above  the  mean  water  level;  and 
z,  the  instantaneous  depth  of  cylinder  immersion.   The  first  and  second 
derivatives  of  n,  with  respect  to  the  time  are  denoted  by  f\   and  r\.      The 
added  mass  is  given  by  Taylor  [_5~]   as: 

I  2r27T      ,1-COS     0N  7T.   .  -..  /       •  «         n\~l  /ox 

m  =  j  Pr  L~3~  (2tt  -  e  }   T(  '  ~cos  9)  +  (sin6-9)]        (2) 

in  which  r  represents  the  radius  of  the  cylinder,  and  9  is  defined  as 
shown  in  Fig.  I . 
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Fig.  I.   Geometry  for  theoretical  analysis  of  wave  impact, 
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The  motion  of  the  free  surface,  r\,    is  related  to  the  maximum  ampli- 
tude by: 

n  =  A  sin  2-rrt/T  (3) 

where  A  and  T  represent  the  amplitude  and  period  of  the  free  surface 

oscillations.   Equation  (I)  can  also  be  written  in  the  form  of  a  slamming 

2 
coefficient,  C  =  2F/pDL  U   as: 
'   s  m 

7  9>"  ,-   7  ,  r   .   2ir  ,  L  3m    2  2tt  ... 

C  =  A .  — rr   -  (m  +  A.)  -r-  s i n  -=-  t  +  —  cos  -=-  t  ( 4 ) 

s    i  , ,  2         i   A     T      ~-      T 
U  dz 

m 

where 

2  2 

A.  =  A./r  ,  m  =  m/pr 
II 

z  =  z/r  ,  U  =  2-itA/T 
m 

—  8 

With  z  =  r(l-cos-»0,  Equation  (2)  can  be  used  to  relate  the  added 

mass  to  the  depth  of  immersion.   Figure  2  is  a  plot  of  Sm/3z  as  a  func- 
tion of  z/D.   Clearly,  9m/8z  begins  with  an  initial  value  of  u  and  drops 

2 
rapidly.   The  quantity  gr/U   is  related  to  the  Froude  number,  Fr,  by 

Fr  =  U  2/2gr.   Thus  C  =  f(Fr,  I  ,  — >  • 

3z 
The  rate  of  change  of  the  normalized  added  mass  with  z  depends  on  0 

and  hence  on  the  time  measured  from  the  instant  of  impact  (See  Fig.  2). 

For  example,  for  the  case  where  n,  =  0  and  t  =  0,  8m/3z  =  tt  and  C  =  tt. 

Consequently,  the  for  the  particular  case  under  consideration,  C  at  the 

instant  of  impact  does  not  depend  on  either  the  size  of  the  cylinder  or 

the  flow  parameters. 
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Fig.    2.      3m/3z   as  a    function   of    z/D. 


16 


The  fact  that  the  cylindrical  members  of  a  structure  in  the  splash 
zone  are  not  necessarily  located  just  at  the  mean  water  level  requires 
the  determination  of  the  particular  value  of  n  for  which  the  slamming 
force  at  the  time  of  impact  is  a  maximum.   It  can  be  demonstrated  through 
the  use  of  Equation  (4)  that  the  maximum  impact  force  occurs  for  the 
case  of  n  =0.   For  this  purpose  Equation  (4)  was  evaluated  with  the  aid 
of  a  computer  program  which  allowed  variations  in  r/A  and  r\   /A.  C  was 
then  plotted  as  a  function  of  z/D  from  zero  to  unity.   Figure  3  is  an 
example  of  one  of  these  plots  and  is  typical  for  most  cases.   The  variables 
were: 

r/A  =  0. 125 

n  /A  =  0.0,  0.4,  and  0.8 

As  can  be  seen  readily  from  Fig.  3,  C   is  largest  at  z/D  =  0  for  n,  /A  =0.0 
and  starts  at  a  value  of  tt  and  drops  to  a  minimum  at  z/D  of  approximately 
0.5.      It  rises  again  as  z/D  approaches  1.0.   Thus  Equation  (4)  indi- 
cates that  C  ,  and  consequently  the  impact  force,  is  of  an  impulsive 
nature  beginning  with  a  finite  value  (not  zero)  at  the  instant  of  impact. 
Since  viscous  effects  are  neglected,  one  would  expect  the  solution  to 
deviate  from  the  actual  situation  as  the  cylinder  becomes  more  fully 
immersed.  Where  this  becomes  the  case  can  only  be  determined  by  experiment. 

However,  one  can  assume  that  the  solution  would  be  valid  for  small 
depths  of  immersion,  which  is  the  case  for  the  instant  of  impact  up  to 

z/D  of  0.02  or  so.   With  this  restriction,  A.  and  m  are  small,  as  is 

2 
sin  2irt/T.  Additionally,  cos  2iTt/T  is  very  nearly  equal  to  1.0.   There- 


fore C  reduces  to; 
s 


C  ~  C°  =  Om/3z),  n 

s    s         t=0 
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Fig.  3.   Theoretical  values  of  C.  as  a  function  of 
z/D  for  several  values  of  n  /A. 
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It  is  not  realistic  to  assume  that  the  impact  force  rises  from  zero 
to  tt  instantaneously.   Several  factors,  specifically  the  compressibility 
of  the  air  between  the  cylinder  and  water  surface,  entrapped  gases  in 
the  water,  and  surface  irregularities  would  account  for  some  finite  rise 
time.   Nonetheless,  the  rise  time  can  be  expected  to  be  short,  i.e.  in 
the  order  of  milliseconds.  The  exact  nature  of  the  rise  is  an  interest- 
ing question  for  further  study.  However,  in  this  analysis  C   is  assumed 
to  vary  I  inearly  during  a  rise  period,  t  .   For  t  >  t  ,  C  drops  from  it 
as  3m/3z.   Figure  4  is  a  representation  of  this  assumption.   Exactly  how 

lonq  of  a  time  interval  t   is  will  be  discussed  later. 
3  r 

The  realization  that  the  impact  force  is  of  an  impulsive  nature  requires 
consideration  of  the  fact  that  this  force  does  not  act  on  a  perfectly  rigid 
body,  but  rather  on  a  cy I i nder  which  is  supported  elastically.   The  re- 
sponse of  such  a  system  approaches  that  of  a  rigid  body  only  if  its 
natural  frequency  approaches  infinity.  Additionally,  the  response  of  the 
system  to  an  impulsive  force  is  heavily  dependent  on  the  exact  nature  of 
the  force  itself  as  well  as  on  the  system  natural  frequency. 

In  general,  the  instantaneous  d isp I acement, x(t), of  a  system  of  mass, 
M,,with  a  spring  constant,  k,  i  s  given  by  (See  e.g.  E6J): 


t 

F(i 
o 


/ 


in  which  £  represents  a  dummy  variable;  F(t),  the  driving  force;  and 

g(t),  the  response  to  a  unit  step  excitation.   Equation  (6)  is  readily 

recognized  as  the  DuHamel  superposition  integral,  which  can  be  expanded 

as  fol  lows  for  q(t)  =  r. —  sin  w  t: 
a      Moo       n 
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Fig.    4.      C    (t)    as  a    function   of   t. 
3  s 
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,t 

"  M^T  C0S  a)nt  /  F(  +  )  sin  aJnt' 
n 

o 


1 

/ 


n        J 


F(t)  cos  a)  tdt 

n 


(7) 


in  which  con  is  the  natural  frequency  of  the  cylinder  and  supports. 

Letting  F(t)  =  F  f(t)  and  changing  the  variable  of  integration  to 


9  =  a)  t,  one  obtai  ns 
n 


k  x(t) 


9/u) 


e/oo 


~^-  =  -  cos  9  /  f  (t)  sin  9d9  +  sin  9  / 


f(t)cos  9d9 


(8) 


in  which  k  =  oon  M,  the  spring  constant.   The  term  on  the  left-hand  side 

can  be  interpreted  as  the  ratio  of  the  instantaneous  force  sensed  by  the 

supports  of  the  cylinder  to  the  actual  mean  force  acting  on  the  cylinder, 

In  other  words,  the  response  of  the  system  at  time  t  is  k  x(t)  which  may 

be  different  from  F  .   For  the  simDle  case  In  which  F(t)  is  a  step 

o  '  r 

function  as  shown  below, 


equation  (8)  reduces  to 


k  x(t) 


=  (I 


-  COS  CO  t) 

n 


(9) 
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It  is  apparent  from  Equation  (9)  that  the  instantaneous  force  sensed  by 
the  system  can  be  anywhere  from  zero  to  two  times  the  actual  mean  force, 
If  one  assumes,  as  before,  that  the  impact  force  is  as  described  by 
Fig.  4,  then  Equation  (8)  must  be  evaluated  using  C  (t)  as  F(t).   Addi- 
tionally, damping  can  be  taken  into  account  by  writing  g(t)  as  L6~]: 

g(t)  =  e'^n1"  sin  oo  t  (  10) 

J  n 

where  C  is  the  damping  coefficient.   Thus  Equation  (8)  can  be  rewritten 

by  replacing  the  forcing  function,  F(t),  by  C  (9)  as: 


>/go  e/w 

n  n 

C  ° 
s 


/  Cs(9)e"^9sin6de  +  cos9  /  Cs(9)e"C9cos9d9 


Equation  (II)  was  solved  with  the  aid  of  a  simple  computer  program  which 

used  a  trapezoidal  integration  of  variable  step.  Values  of  w  were 

taken  as  358  sec   and  628  sec   which  corresponded  to  the  measured  values 

of  co  for  a  6.0  inch  and  a  3.0  inch  diameter  cylinder  respectively,  as 
n 

discussed  later.  The  damping  coefficient  was  0.0I4  in  both  cases,  also 
corresponding  to  measured  values.   The  rise  time  was  varied  from  zero  to 

approximately  0.025  seconds.   Figures  5  through  8  are  representative  plots 

k  x( t ) 
of  for  the  6  inch  cylinder,  for  various  values  of  t  . 

C  ° 

S 
Figure  5  represents  Equation  (II)  plotted  for  a  rise  time  of  0.000I 

second;   Fig.  6,  for  a  rise  time  of  0.0I00  second;   Fig.  7,  for  a  rise 

time  of  0.0I95  second;  and  Fig.  8,  for  a  rise  time  of  0.0230  second. 

k  x( t ) 
Note  that  for  a  very  short  rise  time  the  value  of  reaches  a  value 

C  ° 

5 

of  approximately  I. 7  at  t  =  0.00I  seconds  and  then  drops  off.   As  the  rise 
time  increases,  the  value  of  the  first  peak  drops  off,  with  a  "double 
peak"  appearing  at  a  rise  time  of  0.0I95  seconds. 
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The  interpretation  of  the  above  results  indicates  that,  depending  on 
the  rise  time  of  C  (t),  values  of  the  apparent  slamming  coefficient,  C  , 
may  lie  between  roughly  0.5  and  1.5  of  the  theoretical  instantaneous 
value  of  tt.  Again  this  applies  only  to  a  small  depth  of  immersion,  cor- 
responding to  the  initial  moments  of  impact.   The  significance  of  this 
fact  is  that  values  of  C  measured  in  the  laboratory  may  show  wide  scat- 
ter depending  on  the  rise  time  for  any  given  experiment.   Additionally, 
if  the  surface  is  not  perfectly  plane,  rise  times  may  vary  from  experiment 
to  experiment  resulting  in  an  apparent  non-repeatability. 

Figure  9  is  a  plot  of  k  x(t)/C  as  a  function  of  f  t  for  both  the 
a  r  s  n   r 

first  and  the  second  peak,  usinq  uj  values  for  both  a  6.0  inch  and  a  3.0 

inch  cylinder.   From  this  figure,  the  drop  in  k  x(t)/C  for  the  first 

peak  with  increasing  rise  time  is  apparent.   It  is  also  evident  that  the 

effect  of  ca   is  confined  to  a  narrow  ranqe  of  f  t  values  smaller  than 
n  a      n   r 

about  0.5.   The  force  acting  on  the  cylinder  is  amplified  by  the  dynamic 

response  of  the  system  for  a  value  of  f  t  smaller  than  0.9  and 
r  '  n   r 

attenuated  for  f  t  values  larqer  than  0.9.   The  second  peak  occurs 
n  r  a  K 

only  for  f  t  values  larqer  than  1.09.   The  value  of  k  x(t)/C   increases 
'  n   r  a  s 

from  0.84  to  about  unity  as  f  t   increases  from  1.09  to  1.5. 

'     n   r 
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III.   EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

The  equipment  used  to  generate  the  harmonically  oscillating  flow  has 
been  used  extensively  at  this  facility  over  the  past  two  years.   The 
apparatus  is  described  in  reference  [7]].  The  salient  features, 
as  well  as  the  adaptation  for  this  work,  are  briefly  described  in  the 
fol lowi  ng. 

A.  U  TUNNEL 

Figure  10  is  a  schematic  diagram  of  the  U  tunnel.   A  butterfly  valve 
arrangement  at  the  top  of  the  left  leg  allows  that  side  to  be  completely 
sealed  off.   Compressed  air  can  then  be  admitted  until  the  water  level 
is  driven  to  the  desired  height  in  the  opposite  leg.   Rapid  opening  of 
the  butterfly  valve  allows  the  water  to  oscillate  at  a  natural  frequency 
determined  by  the  geometry  of  the  tunnel.   The  tunnel  contains  approxi- 
mately 5000  gallons  of  water  and  the  period,  T,  is  5.5  seconds.   Ampli- 
tudes up  to  4.0  ft  can  be  generated. 

Previous  work  with  this  tunnel  has  demonstrated  its  capability  in 
studies  involving  oscillating  flow.   Measurements  have  confirmed  that 
separation  does  not  occur  at  the  corners  and  there  are  no  discontinuities 
where  the  individual  sections  are  joined,  which  might  otherwise  lead  to 
disturbances  in  the  flow.  These  factors  were  important  to  this  project 
since  uniformity  of  the  free  surface  was  a  major  concern. 

B.  TEST  CYLINDERS 

Three,  five,  six,  and  six  and  one  half  inch  diameter  aluminum 
cylinders  were  used  in  measuring  impact  forces.   Measuring  2.98  ft  in 
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length,  they  were  constructed  of  aluminum  pipe,  turned  to  the  final 
diameter,  and  polished  to  assure  a  hydrodynamica I  I y  smooth  surface. 
Caps  fitted  to  each  end  prevented  water  leakage  into  the  cylinder. 
Each  cap  contained  a  double  precision  ball-bearing  mounted  flush  with 
the  face.  As  wi  I  I  be  noted  later,  the  force  transducers  were  attached 
to  the  cylinder  via  these  bearings, and  the  purpose  was  to  allow  free 
rotation  of  the  cylinder. 

C.   FORCE  MEASUREMENTS 

Two  canti  lever-beam  force  transducers  were  used  to  measure  the  in- 
stantaneous in-line  (drag)  and  transverse  (lift)  forces.   The  output 
consisted  of  an  electrical  signal  from  strain  gages  attached  to  the 
beam.  The  rounded  end  of  the  transducer  assembly  was  designed  to  fit 
snugly  in  the  bearings  in  the  ends  of  the  cylinders.   Neoprene  "0" 
rings  were  added  to  the  end  of  the  transducers  to  provide  axial  position- 
ing of  the  cylinder  within  the  tunnel. 

Calibration  of  the  transducers  was  accomplished  by  hanging  a  load 
from  the  center  of  the  cylinder.   This  not  only  established  the  level 
of  the  electrical  output  for  a  known  load   but  also  reaffirmed  that  the 
transducers  were  correctly  oriented. 

As  with  other  elements  of  the  experimental  apparatus,  these  trans- 
ducers had  been  used  frequently  at  this  facility  for  over  two  years  and 
the  description  of  their   construction  may  be  found  in  reference  C7U . 
No  changes  have  been  noted  in  the  calibration  of  these  transducers  since 
their  installation,  a  fact  which  establishes  some  degree  of  confidence 
in  their  use. 

Initial  experimental  efforts  established  that  the  transverse  forces 
were  very  small  compared  to  the  in-line  forces, and  consequently  measure- 
ment of  these  forces  was  discontinued  early  in  the  study. 
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D.   MEASUREMENT  OF  FLUID  VELOCITY 

Two  pressure  taps,  one  in  each  leg  of  the  tunnel,  were  connected  by 
tubing  to  a  differential  pressure  transducer  located  midway  between  the 
taps.   The  acceleration  of  the  fluid  can  be  determined  from  Ap  =  ps  -rr  > 
in  which  Ap  represents  the  pressure  differential;  p,  the  fluid  density; 
s,  the  distance  between  the  taps;  and  -rr   ,  the  instantaneous  acceleration 
of  the  fluid.  The  maximum  velocity  and  amplitude  are  related  to  the 
acceleration  by: 


(dU)   =^Lu  =  (2JL)2  a  (12) 

dt  m       m    T 


E.   RECORDING  EQUIPMENT 

A  three  channel  strip  chart  recorder  was  used  to  record  simultaneously 
the  output  of  the  in-line  force  transducer  and  the  differential  pressure 
transducer.  The  third  channel  of  the  recorder  was  also  used  to  record 
the  force  measurements,  but  the  signal  was  first  passed  through  a  low- 
pass  filter  to  remove  the  oscillation  of  the  cylinder  at  its  natural  fre- 
quency, thus  providing  an  "average"  force  record. 

Inasmuch  as  the  impact  force  to  be  measured  consisted  of  a  rather 
short  rise  time,  the  bandwidth  of  the  recorder  was  of  some  importance. 
The  frequency  response  of  this  recorder  was  flat  to  100  Hz.   In  order  to 
verify  the  response  under  the  test  conditions,  several  experimental  runs 
were  made  in  which  the  force  was  both  recorded  using  the  strip  chart  re- 
corder and  a  tape  recorder  simultaneously.   The  taped  signal  was  then 
fed  back  to  the  strip  chart  recorder  at  1/4  speed.   This  method  assured 
that  the  frequency  response  of  the  strip  chart  recorder  was  adequate. 
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F.   EXPERIMENTAL  PROCEDURE 

The  success  or  failure  of  any  experimental  work  depends  on  careful 
attention  to  those  factors  which  may  have  a  significant  effect  on  the 
measurements  being  taken.   For  this  reason, extreme  care  was  taken  with 
all  flow  parameters. 

Prior  to  filling  the  tunnel,  the  cylinder  was  struck  lightly, and  the 

natural  frequency  and  damping  coefficient  were  determined.   The  tunnel 

was  then  filled,  immersing  the  cylinder  several  diameters,  and  again  the 

natural  frequency  was  recorded.   Thus  the  value  of  to  was  obtained  both 

in  air  and  in  water.   In  the  calculations  the  value  of  to   in  air  was 

n 

used  since  at  impact  and  shortly  thereafter,  the  cylinder  behaves  as 
though  it  were  not  immersed. 

Since  the  height  of  the  mean  water  level  with  respect  to  the  bottom 
of  the  cylinder  was  shown  by  the  analysis  to  be  important,  the  water 
level  in  the  tunnel  was  adjusted  carefully  by  slowly  filling  or  emptying 
the  tunnel  until  a  slight  ripple  occurred  at  the  water  surface  due  to 
contact  with  the  cylinder.   This  method  assured  that  the  mean  water  level 
coincided  with  the  cylinder  bottom  surface. 

Of  paramount  importance  was  the  condition  of  the  free  surface  at 
the  instant  of  impact  with  the  cylinder.   Ideally  the  surface  should  be 
perfectly  plane  with  no  surface  disturbances.   It  is  this  requirement 
that  lends  a  certain  degree  of  difficulty  in  the  experimental  procedure. 
Until  the  fluid  in  the  tunnel  has  completed  one  cycle  or  so,  the  period 
of  oscillation  is  not  established.  This  requires  that  the  cylinder  be 
immersed  at  least  once  before  the  period  of  oscillation  has  settled  down, 
resulting  in  some  surface  disturbance.   The  observations  of  the  wave 
impact  indicated  that  one  slam  generated  minimal  surface  disturbances, 
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with  the  disturbances  increasing  with  subsequent  impacts.   For  this 
reason,  all  data  were  recorded  for  the  first  impact  only, and  successive 
impacts  were  discounted.   Approximately  ten  minutes  were  allowed  between 
each  experimental  run  in  order  to  allow  the  water  surface  to  reach  a 
quiescent  state. 

The  stripchart  recorder  was  run  at  a  speed  of  200  mm  per  second,  each 
division  representing  0.005  second.   Figure  II  is  a  typical  example  of 
the  data  obtained.  The  top  trace  is  the  filtered  force  trace, and  the 
bottom  one  is  the  unfiltered  force  trace.   Figures  12  and  13  show  addi- 
tional unfiltered  force  traces  for  various  initial  flow  conditions. 

G.   REDUCTION  OF  DATA 

Two  force  coefficients  were  defined  in  this  investigation.   They  are: 


2F 
C  =  — =—j  (13) 

S   pDLU 
K   m 


and 


C  =-2F   _gD7T  (14) 

pDLU     2U 
m      m 


The  first  of  these  refers  to  the  slamming  coefficient  as  defined  by 
Equation  (4).  The  second  coefficient  represents  the  second  maximum  of 
the  normalized  force  with  the  buoyancy  subtracted.   This  maximum  occurs 
after  the  cylinder  is  fully  immersed. 

As  mentioned  previously,  the  stripchart  recorder  was  run  at  200  mm 
per  second,  thus  each  mm  corresponded  to  0.005  second a   The  measured 
force  was  then  read  each  0,005  second  up  to  0.5  seconds  or  so,  depending 
on  the  diameter  of  the  cylinder.  These  points  were  then  evaluated  with 
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the  aid  of  a  computer  program  which  solved  Equation  (4)  and  plotted 

C  as  a  function  of  z/D.   Figure  14  is  typical  of  such  plots,  which  shows 

the  measured  va I ue  of  C  ,  the  theoretical  value  of  C  from  Equation  (4), 

s  s       M         ' 

and  the  normalized  buoyant  force. 

The  filtered  force  trace  obscured  the  magnitude  of  the  initial  impact, 
However,  after  the  initial  peak,  the  filtered  trace  was  found  to  corre- 
spond to  the  average  value  of  the  unfiltered  trace.  Therefore,  the  data 
points  during  the  early  stage  of  impact  were  read  from  the  unfiltered 
trace, and  the  remainder  were  read  from  the  filtered  trace.   The  main  rea- 
son for  this  is  the  obvious  difficulty  in  interpreting  the  average  value 
which  is  obscured  by  the  oscillation  of  the  cylinder. 
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IV.   DISCUSSION  OF  RESULTS 

The  initial  value  of  the  slamming  coefficient,  C  ,  was  plotted  as  a 

S         Um2 
function  of  various  parameters  such  as  the  Froude  number,  Fr  =  — =—  :  the 

'     go 

Reynolds  number,  Re  =  U  D/v  ;  and  the  Keu legan-Carpenter  number,  K  = 

U  T/D.  These  plots  have  revealed  that  there  is  no  identifiable  corre- 
m 

lation  between  C  and  the  said  parameters.   In  fact,  it  is  the  realiza- 
s  ' 

tion  of  this  lack  of  correlation  that  led  to  the  consideration  of  the  rise 
time,  t  .   Subsequently,  the  amplification  or  attenuation  of  C  was  calcu- 
lated by  dividing  C  by  its  theoretical  value  of  tt.   Then  the  correspond- 
ing f  t  values  were  determined  from  Fig.  9.   Figure  15  is  a  plot  of 

2 
f  t  as  a  function  of  the  inverse  Froude  number  gD/U   .   Even  though 
n   r  3        m  3 

there  is  considerable  scatter,  this  figure  shows  that  for  small  values 
of  U  ,  or  large  values  of  D,  the  rise  time  is  larger  and  thus  yields 
smaller  amplification   (See  Fig.  9.).  For  larger  impact  velocities  or 
for  smaller  cylinders,  the  rise  time  is  considerably  smaller,  and  the 
amplification  factor  approaches  1.7.   |f  should  be  noted 
once  more  that  surface  irregularities,  dissolved  air,  and  other  experi- 
mental uncertainties  do  not  permit  a  more  accurate  determination  of  the 
rise  time.  One  must  add  that  the  rise  time  may  vary  also  with  the 
physical  characteristics  of  the  fluid  above  the  free  surface,  the  angle 
of  inclination  of  the  cylinder,  the  length  to  diameter  ratio,  etc.  Be 
that  as  it  may,  the  Froude  number  appears  to  be  the  dominant  parameter 
as  in  all  free  surface  flow  problems. 

The  physical  reasoning  behind  the  correlation  of  the  normalized  rise 

2 

time, f      t  ,with   qD/U        is  as    follows.      The  disturbances    in   the   free 
n      r'  a       m 
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surface  during  the  first  passage  of  the  flow  increase  with  increasing 
diameter.   On  the  other  hand,  the  larger  the  amplitude,  the  longer  the 
time  for  the  disturbances  to  attenuate.   Thus,  the  rise  time  decreases 

with  increasing  amplitude.   Finally,  the  larger  the  acceleration  of  the 

2 
flow,  i.e.  A(2tt/T)   relative  to  the  gravitational  acceleration,  the 

faster  is  the  rate  of  attenuation  of  the  disturbances  C8,9H.  A  simple 

dimensional  analysis  shows,  therefore,  that  f  t   increases  with 
'  n   r 

2 
gD/A(T  /A).   Since  U   is  proportional  to  A/T,  one  finds  that  f  t   is 

2 
a  function  of  gD/U 
3   m 

A  careful  comparison  of  the  measured  and  calculated  force  traces 
during  the  early  stages  of  impact  lends  ample  credence  to  the  idea  of 
rise  time.   Figures  5  through  8  and  II  through  13  show  that  for  very 
short  rise  times  there  is  a  single  sharp  peak  in  the  force  trace.   For 
intermediate  rise  times  (See  Figures  7  and  12.)  a  double  peak  occurs. 
Finally,  for  larger  rise  times  (See  Figures  8  and  13.)  the  first  peak 
practically  disappears,  and  the  response  to  the  slamming  force  is 
attenuated  by  the  system. 

It  is  event  from  the  foregoing  that  the  determination  of  the  initial 
slamming  force  is  extremely  difficult  because  of  its  magnification  or 
attenuation  by  the  dynamic  response  of  the  system.   The  evidence  presented 
above  shows  that  the  normalized  slamming  force  should  be  assumed  equal  to 
it  and  that  the  response  of  the  system  be  determined  through  a  straight- 
forward vibration  analysis.   Such  an  analysis  will  yield  the  largest 

force  amplification  for  a  normalized  rise  time  of  f  t  —  0.25.   For 
r  n   r 

design  purposes,  the  maximum  slamming  coefficient  may  be  as  large  as  l.7ir. 

Following  the  initial  impact,  the  net  force  acting  on  the  cylinder 
begins  to  decrease  since  the  8m/8z  contribution  decreases.   During 
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this  stage  of  the  flow  the  cylinder  undergoes  damped  oscillations  at  its 
natural  frequency.   As  the  free  surface  rises  further,  the  buoyant  force 
i ncreases,and  the  separation  effects  begin  to  give  rise  to  larger  drag 
forces.  The  buoyancy  subtracted  fluid  force  reaches  its  maximum  at  z/D 
values  from  about  1.5  to  2.0.   Even  though  the  present  flow  situation, 
in  which  there  is  a  free  surface, cannot  be  directly  compared  with  an 
impulsive  flow  about  a  circular  cylinder,  the  rise  of  the  drag  force  to 
a  maximum  at  z/D  —  1.75  is  very  much  like  the  rise  of  the  drag  coeffi- 
cient to  a  maximum  at  a  relative  fluid  displacement  of  2.0  in  impulsive 
flow  [J  0]].  The  said  rise  in  the  drag  coefficient  is  because  of  the  for- 
mation of  two  symmetrical  vortices  behind  the  cylinder.  As  the  motion 
continues, the  vortices  become  asymmetrical  and  shed  a  I ternat i ng I y. 

The  maximum  force  corresponding  to  the  second  peak  has  been 

I    2 
corrected  for  buoyancy,    normalized  by  -^   pU  DL,and  plotted  as  a 

function  of  the  Froude  number  in  Figure  16.   Two  facts  are  apparent: 
First,  the  scatter  in  the  data  is  considerably  less  than  that  correspond- 
ing to  the  initial  impact;  second,   C,  [See  Equation  (14).]  decreases 
with  increasing  Froude  number  to  a  value  of  about  unity  and  remains 
nearly  constant  for  all  Froude  numbers  from  0.5  to  2.0. 
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V.   CONCLUSIONS 

The  theoretical  and  experimental  investigation  of  flow  impact  on 
circular  cylinders  warranted  the  following  conclusions: 

1.  The  dynamic  response  of  the  system  is  as  important  as  the  impact 
force  and  one  cannot  be  determined  without  taking  the  other  into 
consi  deration . 

2.  The  initial  value  of  the  slamming  coefficient  is  essentially  equa 
to  its  theoretical  value  of  it,  and  the  system  response  may  be 
amplified  or  attenuated  depending  upon  its  dynamic  characteristic, 

3.  The  normalized  rise  time  appears  to  be  a  function  of  the  Froude 
number. 

4.  Following  the  impact,  the  cylinder  undergoes  damped  oscillations 
at  its  natural  frequency. 

5.  The  buoyancy  corrected  normalized  force  in  the  drag  dominated 
region  reaches  a  maximum  at  a  relative  fluid  displacement  of 
about  z/D  —  1.75.   Subsequently,  the  shedding  of  vortices  as 
well  as  the  deceleration  of  flow  reduces  the  maximum  drag  coeffi- 
cient to  below  unity. 

6.  Since  it  does  not  appear  that  the  rise  time  can  be  determined 
theoretically,  the  dynamic  response  of  the  system  should  be 
analyzed  using  the  theoretical  value  of  the  slamming  coefficient 
and  the  experimentally  determined  rise  times. 
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